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^ ■ Complex potentials are constructed as Darboux-deformations of short range, 

[t I . radial nonsingular potentials. They behave as optical devices which both refracts 

I and absorbs light waves. The deformation preserves the initial spectrum of energies 

' and it is implemented by means of a Gamow-Siegert function (resonance state). As 

■ straightforward example, the method is applied to the radial square well. Analytical 
i derivations of the involved resonances show that they are 'quantized' while the 

^p^l corresponding wave-functions are shown to behave as bounded states under the 

■ broken of parity symmetry of the related one-dimensional problem. 

: 

1 Introduction 

I ■ 
> • 

CN , Solutions of the Schrodinger equation associated to complex eigenvalues e = and 

^ I satisfying purely outgoing conditions are known as Gamow-Siegert functions [1,2]. These 

Tj- ■ solutions represent a special case of scattering states for which the 'capture' of the incident 

fS) ■ wave produces delays in the scattered wave. The 'time of capture' can be connected with 

^ , the lifetime of a decaying system (resonance state) which is composed by the scatterer 

! and the incident wave. Then, it is usual to take Re(e) as the binding energy of the 

^ I composite while Im(e) corresponds to the inverse of its lifetime. The Gamow-Siegert 

functions are not admissible as physical solutions into the mathematical structure of 
Quantum Mechanics since, in contrast with conventional scattering wave-functions, they 
, are not finite at r — > oo. Thus, such a kind of functions is acceptable in Quantum 

Mechanics only as a convenient model to solve scattering equations. However, because of 
the resonance states relevance, some approaches extend the formalism of quantum theory 
so that they can be defined in a precise form [3-9]. 

In this paper the Gamow-Siegert functions are analyzed not to precisely represent 
decaying systems but to be used as the cornerstone of complex Darboux-transformations 
in Supersymmetric Quantum Mechanics (Susy-QM). The 'unphysical' behaviour of these 
solutions plays a relevant role in the transformation: While the Gamow-Siegert function 
u diverges at r — >■ oo, the limit of its logarithmic derivative (3 = —u'/u at r ^ cxd is a 
complex constant. This last function is used to deform the initial potential into a complex 
function. It is worth noticing that complex eigenvalues e have been used as factorization 
constants in Susy-QM (see for instance [10-18] and the discussion on 'atypical models' 
in [19]). However, as far as we know, until the recent results reported in [20-22] the 



connection between Susy-QM and resonance states has been missing. Our interest in 
the present work is two-fold. First, we want to show that appropriate approximations 
lead to analytical expressions for the real and imaginary parts of the resonance energies 
in the case of radial short range potentials. Second, we want to enlarge the number of 
applications of the Darboux transformations in constructing new exactly solvable models 
in Quantum Mechanics. 

In Section 2 the main aspects of solving the Schrodinger equation for radial, nonsin- 
gular short range potentials are reviewed. The relevance of the scattering amplitude S{k) 
in the analysis of both the bounded and scattering wave-functions is clearly stated. In 
Section 3 the analytical properties of S{k) as a function of the complex kinetic parameter 
k are studied. It is shown that poles of the scattering amplitude which live in the lower 
fc-plane are connected with the Gamow-Siegert functions while the poles on the positive 
imaginary axis lead to bounded physical energies. The general aspects of the complex 
Darboux-transformations are included in Section 4. All the previous developments are 
then applied to the square radial well and new analytical expressions for the involved res- 
onance energies are reported in Section 5. We conclude the paper with some concluding 
remarks in Section 6. 



2 Short range radial potentials revisited 



2.1 General considerations 



Let us consider the Hamiltonian of one particle in the external, spherically symmetric 
field U{r) = U{r). The Schrodinger equation reduces to the eigenvalue problem: 
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A + f/(r) ^(f) = S^lj{r). 



We are interested in potentials which decrease more rapidly than 1/r at large r. Indeed, 
U{r) is a short-range potential (i.e., there exists a G Dom(f/) := Djj such that U{r) = at 
r > a) which is nonsingular (i.e., U{r) satisfies r'^U{r) — at r ^ 0). The wave function 
%l){r) must be single- valued and continuous everywhere in Du. Since the operators if, L^, 
L3 and P commute with each other, in spherical coordinates the Schrodinger equation ([T]) 
has solutions of the form 

V'£m(r) =i?^(r)y,„(^,(^) (2) 

where 9 and are respectively the polar and the azimuthal angles, R(,{r) is a function 
depending on r and Yf^iO^ip) stands for the spherical harmonics. For simplicity in cal- 
culations we shall use the function ui{r) = rR^{r). The introduction of ([2]) into equation 
yields: 
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i{i + 1) 



+ f (r) 



ue(r 



+ VAr) 



ue{r) = k'^UfXr) 



(3) 



with V = 21711/ /h'^ and k"^ = E = 2mEIY? the dimensionless expressions for the potential 
and the kinetic parameter respectively. As usual, the effective potential Vi is defined as 
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the potential v plus the centrifugal barrier + l)/r^. Physical solutions ip{j) of ([T]) 
should be constructed with functions U(, fulfilling the following boundary conditions: 

(4) 

Furthermore, the function and its derivative have to be continuous for r > since 
([3]) includes second order derivatives. As we shall see, the study of the involved matching 
conditions is simplified by introducing a function /?£ as follows 

A(r):=-^ = -;^lnu,(r). (5) 

Ui[r) dr 

In general we will assume that the energy spectrum is composed by negative (bounded) 
and positive (scattering) energies. For the sake of notation we shall use the symbol 0£(r) 
for the solutions of ([3]) which satisfy the boundary conditions (jlj); the symbol U(,{r) will 
stand for general mathematical solutions of equation ([3]). 



2.2 Bases of solutions 

For arbitrary £ and r > a, the appropriate basis of solutions can be written in the form 

u}^\r) = ikrh^l\kr), u\ \r) = —ikr hf'\kr) (6) 

where /li^'' and h^^^ are the spherical Hankel functions of order u = i (see e.g. [23]). For 
large values of kr in the region r > a, the functions ([H]) behave as follows 

(7) 



uf^\r) ^ exp 



±i i kr — — 
2 



Thereby, u^^^ represents a diverging wave and describes particles moving with speed d = k 
in all directions from the origin. In turn, the function u\ represents a converging wave 
and describes particles moving with speed d towards the origin. 

It is worth noticing that the Schrodinger equation ([3]) is invariant under the change 
k —k. Hence, if UfXr, k) is solution of ^ ioi E = k'^, then ug^r, —k) is also a solution 
for the same energy. However, these last functions must differ at most in a constant factor 
because the solutions of ^ have to be single- valued. In contrast, up to a global phase, 
they interchange their roles at r ^ oo, as indicated by the relationship 

uf\r,k) = e'^^uf\r,-k), r ^ oo. (8) 

On the other hand, the basis of solutions at r < a may be written 



uf\r) = ^ krjfXkr), uf\r) = -Ae kr ng{kr), Ai := ^^^2) 

where je and ni are respectively the spherical Bessel and Neumann functions of order i. 
For very small values of r we have 

u^^\r ^ 0) ^ {krY+\ uf\r ^ 0) ^ (10) 
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2.2.1 Bounded states. 



Imaginary values of the kinetic parameter k = ±i^y\E\ = ±iK correspond to negative 
values of the energy E = k"^ = —k^. As usual, we shall consider k to be in the upper 
complex /c-plane /+. In this way u\ \r, k) does not satisfy the boundary condition at 

(2) 

r — >■ cx) and (r, k) does not fulfill the boundary condition at r = 0. Hence, the physical 
solutions should be constructed with uf\r, k) and u^i^\r, n). We write: 

uf^ (r, k) gi,{r, k) r < R 

(11) 

W^)u^^\r,K) r>R 

where the intermediary function gi{r, k) is equal to 1 at r = and depends on the 
potential. The coefficients and are complex numbers while the matching 

condition f3£{a,K) = k, corresponds to a transcendental equation which is fulfilled for 
a set of discrete values of the kinetic parameter Kn, n = 0,1,2, .. . Under these conditions 
the functions (j)p( finite and the integral 

/• + 00 

/ \(t),{T,Kn)\^dr (12) 

converges. Thereby, |0^(r, /t„)p can be identified as the involved Born's probability den- 
sity. 



2.2.2 Scattering states. 

Positive energies E = k"^ correspond to real values of k and both solutions ([7]) remain finite 
in r > a. As a consequence, both of them are physically acceptable in this region. To 
analyze these solutions first let us consider the case i = and, for the sake of simplicity, 
let us drop this subindex from the functions ue{r). The general solution in the free of 
interaction zone (r > a) may be written as 

Mout(r) = 7(^) [u^-\r) - S{k)u^+\r)] . (13) 

The constant '~f{k) is the amplitude of the converging wave and does not depend on 
the potential v{r). The scattering amplitude S{k), in contrast, strongly depends on 
the potential and encodes the information of the scattering phenomenon. The matching 
condition between Mout('") and the solution Uin{r) in the interaction zone reads 

/3out(a) = An (a) (14) 

and can be satisfied by appropriate values of 7 and S. Indeed, equation (IT^ is equivalent 
to the following matrix array 
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the solution of which is given by the system of equations 
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r=a 



(16) 



2ik 

where W{-, ■) stands for the Wronskian of the involved functions. 

To get a better idea of the roles played by each one of the terms in equation (fT3|) let 
us consider the free motion {v{r) = 0) in the whole of D^. From equations ([6]) and (Q we 
realize that a regular solution at the origin may be written as follows 

^frcc('") = u^~\r) — u^~^\r) = —2ikrjQ{kr) = —2ism{kr). (17) 

This term can be introduced into equation (fT3|) to get 

Mout(r) = l{k) [Mfree(r) - (Sik) - l)MW(r)] . (18) 

Thus, the total wave function after scattering Uout is given by the wavepacket one would 
have if there were no scattering 7Mfree plus a scattered term —7(6' — l)u''^\ This last 
example shows the important role played by the scattering amplitude in the analysis of 
positive energies. 

On the other hand, the potential we are dealing with is neither a sink nor a source of 
particles, then the condition |5'(/i;)| = 1 (elastic scattering) is true and we have 

S{k) = e^'^^'l (19) 

The scattering phase 6{k) represents the difference in phase of the outgoing parts. The 
situation is clear if one introduces ([T7I) and f|T9|) into equation f|T8l) to get 

Mout(r ^ 00) = -2t-f{k) e'^^^^ sin(fcr + 5). (20) 

As we have indicated, the scattering states include k E M. since E > Q. However, the 
function f|T3|) is also valid if k is one of the imaginary values of the kinetic parameter 
leading to the bounded energies En = —i^n- Indeed, this case provides further information 
about the scattering amplitude. For fc„ = iHn, the function f[T5]) ai r —> 00 reads 

Uout{r 00, /t„) 2i7(m„) [e^"'' - 5(m„)e-''"''] . (21) 

The first term in square brackets never vanishes («;„ > 0), then the coefficient S must be 
singular at the complex point /c„ = m„. In other words, the points fc„ are nothing but 
poles of the scattering amplitude. Since these points are in the upper complex fc-plane 
/+, we realize that S{k) is a regular function in /+ except at kn, n = 0, 1, 2, . . . We shall 
use these results in the next sections. 

Finally, similar expressions can be found for arbitrary azimuthal quantum numbers i. 
A straightforward calculation gives 

ue{r) = 7,(A;) (u['\r) - Si{kW^\r)^ , r>a (22) 

the asymptotic behaviour of which reads 



Ui[r 00 J 



-2i7£(fc)e*^^('=) sin (^kr + 6e{k) - . (23) 

The above discussed properties of the scattering amplitude S{k) are easily generalized to 
the case of arbitrary angular momentum Se{k). 
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3 Gamow-Siegert functions 



3.1 Analytic continuation of the scattering amplitude Si{k) 

In the previous section we realized that Se is a regular function in /_(_ except at the points 
k = in leading to bound states of the energy. In order to get a well behaved scattering 
amplitude in the whole of the complex fc-plane, this function must be extended to be 
analytic in /_ . With this aim, first let us construct an arbitrary linear combination of the 
functions ([7]), it reads 



ueir) = Ceik)4-\r) - Uk)u[^\r) = Qik) [u^'ir] - Seik)u',''> (r)^ . (24) 

This function is regular at the origin if ^^(r = 0) = 0, then we arrive at the following 
relationship 

Uk) _ u'f ^r) 



Se{k) 



(25) 



Ci{k) 4+)(r) 

Now, after a change of sign in k, the coefficients of the single-valued function (12^ should 
satisfy 

- e'-'U-k) = aCeik), e'^'OX-k) = -a^k) (26) 

where a is a proportionality factor and we have used equation ([8]). In this way, the 
expression (!25|) leads to 

Se{k) Sei-k) = 1. (27) 

Let ki G /+ be a zero of Se- From equation fl27|) we notice that the scattering amplitude 
is well behaved in /_ except at the point —ki, for which Se{—ki) oo. Thus, —ki G /_ 
is a pole of 5*^. 

On the other hand, if the kinetic parameter k is real, then the complex conjugate of 
any solution of equation ([3]) is also a solution for the same energy. However, the solution 
is unique so these last functions differ at most in a global phase. Thereby, in a similar 
form as in the previous case, we get 

Se{k)Se{k) = l (28) 

where the bar stands for complex conjugation. This last equation is a consequence of the 
elastic scattering we are dealing with since |5'£(/c)| = 1. However, as k is real, also equation 
( !28l) has to be extended to the whole of the complex /c-plane. The natural condition reads 

Se{k)Se(k) = l. (29) 

It is clear that Si{ki) diverges because ki G /+ is a zero of Se- In other words, ki G /_ is 
a pole of Se since \Se{ki)\ oo. Moreover, as —ki is a pole of Se, from equation ( 129|) we 
also notice that Se{—ki) = 0. Thus, —ki is another zero of the scattering amplitude. In 
summary, if ka (ki) is a pole (zero) of Se, then —ka {—ki) is another pole (zero) while ka 
and —ka {ki and —ki) are zeros (poles) of the scattering amplitude. In this way Se{k) is 
a meromorphic function of the complex kinetic parameter k, with poles restricted to the 
positive imaginary axis (bound states) and to the lower half-plane /_ (see Figure [I]). 



6 



< 
I 
1 

1 

— ki 

° 


> IK„ 

> 

1 

) 

ka ~ ki 

° 


. • • 

ka = ki 


• • . 



Figure 1: Schematic representation of the poles (disks) and the zeros (circles) of the scattering amplitude 
Si{k) in the complex fc-plane. Bounded energies En — — correspond to the poles located on the positive 
imaginary axis. 



3.2 Resonance states 

Let us express one of the poles of the scattering amphtude as ka = C(i — ia2i with ai an 
arbitrary real number and 0^2 > 0. The function fl22|) . evaluated for /cq, at r — > 00 reads 

u,{r ^ 00, k^) ^ 7,(fc„)e-^('^^^-^-/2) ^^-a,r _ s,{k^)e^^'] . (30) 

Observe that the first term in brackets vanishes for large values of —a2r. Hence, there 
remains only the scattering, purely outgoing wave e"^*^. As a consequence, the particle 
is maximally scattered by the potential field v{r) and the wave function is not "well 
behaved" since it does not fulfill the boundary condition at r ^ 00. However, as we have 
seen, this last "unphysical behaviour" is not only natural but necessary to study the elastic 
scattering process of a particle by a short range, nonsingular radial field interaction. Wave 
functions behaving like ( !30ll are known as Gamow-Siegert solutions of the Schrodinger 
equation ([3]). We write 

{Oi ge{r, ka) uf^ (/c„r) r < R 

(31) 
-li{ka)Se{ka)u['^\r, k^) R<r 

where 61 is a constant and the intermediary function ge{r, ka) is equal to 1 at r = 
and depends on the potential. Besides the matching condition ( fT6l) . the Gamow-Siegert 
functions fulfill the purely outgoing boundary condition: 

lim \ I3f^\r) + tka] = 0. (32) 

r— >+oo L J 

Our interest is now addressed to the Gamow-Siegert functions not precisely as describing 
a decaying system but as appropriate mathematical tools to study the construction of 
complex Darboux-deformed potentials. 
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4 Complex-Darboux deformations of radial poten- 
tials 



4.1 General considerations 

In order to throw further hght on the function /^^(r) we may note that (jSj) transforms the 
Schrodinger equation into a Riccati one 

-(3',{r)+f3^{r) + ee = Ve{r) (33) 

where the energy k'^ has been changed for the arbitrary number e^. Remark that 0331) is 
not invariant under a change in the sign of the function beta: 

Peir) + Plir) + Q = V^ir) + 2/5;(r). (34) 

These last equations define a Darboux transformation Vi{r) = Vi{r,ei) of the initial 
potential Vf(r). If is a nontrivial complex number, then the Darboux-deformation is 
necessarily a complex function 

Ve{r) = 14(r) + 2/3;(r) = r,(r) - 2^ ln(^,(r) (35) 

where the transformation function Uf{r,e) = ipe{r) is the general solution of ([3]) for the 
complex eigenvalue e^. The solutions ye{r) = yeXr, ee, £) of the non-Hermitian Schrodinger 
equation 

-yl{r) + V,{r)y,{T)=£y,{r) (36) 

are easily obtained 

W((^,(r),M^(r)) 
Vfin oc , (37) 

where u^ij-) is eigensolution of equation ([3]) with eigenvalue £. 



4.2 New exactly solvable optical potentials 

Let the transformation function v?e(r) be a Gamow-Siegert function (I3T!) . The superpo- 
tential (3f^^\r, ka) and the new potential Vi{r) behave at the edges of Dy as follows 

I — z/Cq, r — > +00 I Vi[r) r — > +cx3 

It is remarkable that even for an initial potential with i = the corresponding Darboux- 
deformation contains a nontrivial centrifugal term in the interaction zone r < a. Notice 
also that the new potential is mainly a real function at the edges of D^. Thus, the function 
Im(V^) is zero at the origin and vanishes at r ^ +cx3. In general Im(V£) oscillates along 
the whole of (0, +oo) according with the level of excitation of the Gamow-Siegert solution 
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Ui (r, ka) '■ The higher the level of excitation the greater the number of oscillations. Such 
a behaviour implies a series of maxima and minima in Im(V£) which can be analyzed in 
terms of the optical bench [19]. 

As regards the physical solutions of the new potential it is immediate to verify that, 
in all cases (bounded and scattering states of the initial potential), the corresponding 
transformations fulfill the boundary condition at r = 0. The analysis of the behaviour at 
r — >■ oo is as follows. 

4.2.1 Scattering states. 

Let ue{r) be a scattering state of the initial potential in r > a. At r +oo, the 
corresponding Darboux-Deformation ( |371) behaves as follows 

lim ye{r) = [{K + k)u'^-\r) - {K - k)S,{k)u^^\r)\.^ (39) 

where we have used equation (!22ll . As we can see, the Darboux-deformations of behave 
as scattering states at r > a. 

If the scattering state ui{r) is a Gamow-Siegert function and k ^ ka, ka, —ka, —ka, 
the Darboux-deformation preserves the purely outgoing condition (1321) . Special cases are: 

1. k = ka ^ the Wronskian in (1571) is zero and there is no transformation. 

2. k = ka ^ Si{ka) = 0, then yi{r) is an exponentially decreasing function. 

3. k = —ka =^ Se{—ka) = and the coefficient of u^~\r) vanishes, then ?/^(r) = 0. 
A. k = —ka =^ SfX—ka) — * oo, then yt{r-) fulfills the purely outgoing condition (1321) . 

4.2.2 Bounded states. 

If /+ 3 A; = m = i\/\E\ the function uiijr, ^J\E\) corresponds to one of the bounded states 
(fTT]) . Thereby, from equation (1371) we get 

lim y,{r, ^/\E\) ~ {iK - VW\) UVW\) e~^^ e'V^^ (40) 

r— ^+00 

Thus, the Darboux-deformations of bounded wave-functions vanishes at r ^ oo. However, 
although this new functions are square-integrable, they do not form an orthogonal set in 
the Hilbert space spanned by the initial square-integrable wave-functions [14] (see also [24] 
and the 'puzzles' with self orthogonal states [25]). 

In the next section we are going to solve the Schrodinger equation for the radial square 
well. We shall get the solutions corresponding to complex energy eigenvalues and purely 
outgoing boundary conditions. Thus, we shall construct the involved Gamow-Siegert 
functions. 
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5 Gamow-Siegert states of the radial square well 



Let us consider the potential 



f (r) 



—Vq r < a 
a < r 



(41) 



with a and vq positive real numbers. The regular solution of the Schrodinger equation 
for this potential may be written 



•je{qr) 

-i'jeik) kr \hf\kr) + SfXk)hf\kr) 



r < a 
a < r 



(42) 



where the interaction parameter q is defined by = Vq + k"^. For simplicity, we shall 
analyze the s-wave solutions {i = 0) for which the matching conditions (fT6|) lead to 



S{q) 



ik sin(ga) + q cos(ga) 
iksin{qa) — gcos(ga) 



-2ika 



l{k) 



2ik 



[iksm{qa) — q cos (qa)]. (43) 



Let us take 6 = 2ik. The scattering amplitude S{k) has zeros and poles which respectively 
correspond to the roots of the transcendental equations 



= tan(ga), 

k 



~ = tan(ga). 
k 



(44) 



Indeed, equations fl44j) correspond to the even and odd quantization conditions for negative 
energies if k = ik E /+, just as it has been discussed in the previous sections. However, we 
are looking for points k^ in the complex fc-plane such that Iie{ka) 7^ and lm{ka) 7^ 0. 
Thereby, if ka is a pole of S, the converging wave in (l42l) vanishes. Then the Gamow- 
Siegert function reads 



2ikaSm[q{ka)r] r < a 
2ika sm[q{ka)a] e*'^"('"~") a < r 
In order to solve the second equation in flH|) let us rewrite the scattering amplitude as 

-^2iqa(^f^ + g) + (g - fc) 



(45) 



Sik) = - 



j,2iqa(J^ - g) - (g + fc)_ 

In this way, the poles of 5* are also the solutions of the transcendental equation 

2iqa 



\k 


2 _ 




2 + 2Re(fcg) 




k 


- q 


2 



(46) 



(47) 



If the interaction parameter g is real, then sin(2ga) = and we obtain the following 
quantization rule for the energy 



Tin 
2^ 



Er, 



n = 0,l, 



(4J 
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On the other hand, let k and q be the complex numbers k = Ki + iK2 and q = Qi + iQ2- 
Then q"^ = vq + e = (fo + ei) + ie2, with e = k^. Hence we have 



Ql-Ql = vo + ei, 2Q1Q2 = £2 (49) 

and 

Kf - Kl = ei, 2K1K2 = 62. (50) 

The imaginary part of equation fj47|) is sin(2(5ia) = and we obtain the following 
quantization rule 

TTTl 

Qi,n = ^, n = 0,l,... (51) 

Therefore, if \Q2\ « 1, the first equation in (H^ shows that equation leads to 
ei ~ En- The second equation in (H9l) . on the other hand, gives a value of €2 which is not 
necessarily equal to zero. In this way, we shall study the poles of 5* for which the complex 
points e are such that Re(e) ~ En and Im(e) ~ — ^, with < F << 1. The main problem 
is then to approximate the adequate value of F as connected with the discreteness of Qi 
in fl3T]) . With this aim, for a given finite value of a we take \Q2\a. << 1. A straightforward 
calculation shows that the second equation in uncouples into the system 



Qi = ^/v^sm{Qla) cosh{Q2a), Q2 = ^/v^cos{Qla) smh{Q2a). (52) 
In turn, the second one of these last equations reduces to 

cos(Qia) ^ ^— = -. (53) 
ay/vo 7] 



Now, let |(5| << 1 be a correction of Qi around the quantized values (ISTll . that is Qia ^ 
Qi,nO' + 5- Then, equation ( |53l) leads to 

sin(7m/2) 

6n = ^ — —, n odd. (54) 

Thus, only odd values of n are allowed into the approximation we are dealing with. Now, 
to ensure small values of Q2 we propose = X]m=i ^m/ff^- From fH9l) we have 

^ ^ f 0'Qi,n \ ^ _ 2aQi,„ sin(7rn/2) ^ 1 - _ 2A1A2 _ 

vq \ V J V'^ 

In order to cancel the term including r^"^ we take Ai = ±1 (the appropriate sign will be 
fixed below). Higher exponents in the power of l/?7 are dropped by taking Xm>i = 0. In 
order to get ei En, after introducing ( jSTl) into (j55l) . we realize that the condition r] » 1 
allows us to drop the term including ri~^. Thus, at the first order in the approximation 
of Q2a we get ei ~ En- 
It is also necessary to have regard to the positiveness of energy ei, which is ensured 
whenever n exceeds a minimum value. Let us take n := riinf + m, m G Z^, where ninf is 
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the ceiling function of 217/71, i.e. n 



inf 



2a,Jvo 



even (odd) if riinf is odd (even). We finally arrive at 



Since n is odd, the integer m is 



Re(e) 



(rijnf + m)7r 



1 ^^0, 



m 



0,2,4, 
1,3,5, 



ninf odd 
ninf even 



(56) 



On the other hand, it is worth noticing that the introduction of = and fl5^ into 

the second equation in ( l49l) leads to the conclusion that = +I/77 has to be dropped. 
The expression for the imaginary part of the complex energy e is calculated as follows. 
First, let us rewrite equations ( l52i) in the form 

Ki = ±y/vosm{Qia) smh.{Q2a), K2 = ±y/v^cos{Qia) cosh{Q2a). (57) 

Since €2 = — r/2, from equations ( l50ll we realize that the sign of Ki must be the opposite 
one of K2. The introduction of ( l53l) into the second equation in (1571) leads to 



1 

±- 

a 



(5J 



where we have used Q2a ^ —^/v- If V » 1 S^t aK2 ~ ±1. Then, this last result 
together with the real part of the factorization constant e gives 



ei + 



(59) 



Re (k) 




|S(k) I 



Figure 2: The absolute value of the scattering amplitude S{k) close to the pole ka (and the zero ka) 
which corresponds to the first resonance Ea = fc^ of the square well reported in Table [1] for ninf = 64. 
The limits of our approach give the number ka = 2.063412 — i 0.099882. 



Thus, for a^ei > 1 we obtain Ki ^ ±y^. Finally, the introduction of ( !58l) and ( !59l) into 
the second equation of ( l50l) produces 

Im(e) = -^ = -mil \K2\ ~ (60) 

In Figure [2] we show the absolute value of the scattering amplitude close to one of its poles 
ka and the corresponding zero ka- Some of the first resonances are reported in Table [T] 
for different values of the potential strength vq and the cutoff a. 
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5.1 Complex Darboux-deformations of the radial square well 



Figure [3] shows the global behaviour of a typical Gamow- Sieger t function associated with 
the radial square well. Notice the exponential growing of the amplitude for r > a. This 
kind of solutions are used in fl35|) to complex Darboux-deform the square well as it is 
shown in Figure HI 




5 10 15 20 25 5 10 15 20 25 



Figure 3: The real (left) and imaginary (right) parts of the Gamow-Siegert function (|45|) associated 
with the first resonance state reported in Table [T] for Tijnf = 64. 



The cardiod-like behaviour seems to be a profile of the complex potentials derived by 
means of Darboux-deformations (compare with [20] and [22]). In this case, the complex 
potential shows concentric cardiod curves for values of r inside the interaction zone, one 
of them is shown at the left of Figure HI The real and imaginary parts of Vi{r) are then 
characterized by oscillations and changes of sign depending on the position in r < a. As 
a consequence, these complex potentials behave as an optical device which both refracts 
and absorbs light waves (see details in [22] and discussions on the optical bench in [19]). 
The presence of this kind of oscillations is also noted at distances slightly greater than the 
cutoff r = a. Hence, the complex Darboux-deformations ( l35l) are short range potentials 
which enlarge the initial 'interaction zone' as it is shown in the right part of Figure |H 



3000 
2000 
1000 


-1000 
-2000 
-3000 

-5000-4000-3000-2000-1000 -1 -0.5 0.5 1 1.5 

Figure 4: The Argand-Wessel diagram of the complex Darboux-deformed square well with uq = 100 
and a = 10. Left: Detail of the cardiod-like behaviour of the new potential between r = 9.8 (disk) and 
r = 10 (circle). Right: The disk is evaluated at r = 10.1 and the circle at r = 13. As complementary 
information: Va{r = 10) = 100.2255 - i 0.8076 and Vo{r = 10.1) -0.1157 - i 0.8306. 
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ninf = 64 


ninf = 202 


m = 1 
m — 3 

171—5 

m — 7 


04.247696 -i 0.412198 
10.761635 -i 0.656098 
17.472966 -^0.836013 
24.381689 -i 0.987556 


16.791319 -i 0.819544 
36.925312 -i 1.215324 
57.256697 -i 1.513363 
77.785474 -i 1.763921 



Table 1: The first four resonance energies for the radial square well of intensity wo = 100 and cutoff 
a = 10 (ninf = 64) and for the radial square well of intensity = 1000 and cutoff a = 10 (ninf = 202). 
Notice that, in each case, the even values of m obtained for the related one-dimensional problem are 
missing (see e.g. Table A.l of reference [22]). 

6 Concluding remarks 

We have studied the elastic scattering process of a particle by short range, nonsingular 
radial field interactions. The poles of the involved scattering amplitude S{k) play a 
relevant role in the construction of physical solutions connected with either bounded or 
scattering states of the energy. It has been shown that poles ka in the lower half complex 
/c-plane lead to the "unphysical" Gamow-Siegert functions which are necessary to depict 
the behaviour of the scattering phenomenon. 

The Gamow-Siegert functions were used to transform the radial short range potentials 
into complex ones which preserve the initial energy spectrum. These new potentials are 
'opaque' in the sense that they simultaneously emit and absorb flux, just as an optical 
device which both refracts and absorbs light waves. The transformation preserves the 
square-integrability of the solutions at the cost of producing non-orthogonal sets of wave- 
functions. It is also notable that scattering states are transformed into scattering states 
while deformed Gamow-Siegert functions can be either a new Gamow-Siegert function or 
an exponentially decreasing function depending on the involved kinetic parameter k and 
the pole ka- 

The method has been applied to the radial square well in the context of s-waves (£ = 0). 
Analytical expressions were derived for the corresponding complex energies in the long 
lifetime limit (i.e., for small values of |Im(£^Q)| = r/2) and, as a consequence, such energies 
fulfill a 'quantization rule'. In contrast with the 'even' and 'odd' resonances of short range 
potentials defined in the whole of the straight-line (see e.g. [22]), the resonances of radial 
short range potentials are labelled by only odd (positive) integers. This result enforces 
the interpretation of Gamow-Siegert functions as representing quasi-bounded states: The 
bounded spectrum of a potential which is invariant under the action of the parity operator 
v{x) — v{—x) includes odd and even functions. The symmetry is broken by adding an 
impenetrable wall at the negative part of the straight line and only the odd solutions are 
preserved. As we have shown, the same is true for resonance states. 

Finally, the results reported in this paper are complementary to our previous work 
[22]. It is remarkable that explicit derivations of Gamow-Siegert functions are barely 
reported in the literature, not even for simple models like those studied here (however 
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see [5,20,21,26-29]). We hope our approach has shed some hght onto the solving of the 
Schrodinger equation for complex energies and functions fulfilling the purely outgoing 
boundary condition. 
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